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ABSTRACT 

We present the first results of the expected variations of the molecular line emission arising from 
material recently affected by C-shocks (shock precursors). Our parametric model of the structure of 
C-shocks has been coupled with a radiative transfer code to calculate the molecular excitation and 
line profiles of shock tracers such as SiO, and of ion and neutral molecules such as H 13 CO + and 
HN 13 C, as the shock propagates through the unperturbed medium. Our results show that the SiO 
emission arising from the early stage of the magnetic precursor typically has very narrow line profiles 
slightly shifted in velocity with respect to the ambient cloud. This narrow emission is generated in 
the region where the bulk of the ion fluid has already slipped to larger velocities in the precursor as 
observed toward the young L1448-mm outflow. This strongly suggests that the detection of narrow 
SiO emission and of an ion enhancement in young shocks, is produced by the magnetic precursor of 
C-shocks. In addition, our model shows that the different velocity components observed toward this 
outflow can be explained by the coexistence of different shocks at different evolutionary stages, within 
the same beam of the single-dish observations. 

Subject headings: ISM: clouds — ISM: jets and outflows — physical data and processes: shock waves 
- ISM: individual (L1448) 



1. INTRODUCTION 

Shocks are an ubiquitous phenomena in the inter- 
stellar medium. In molecular dark clouds where the 
frac tional ionisat i on o f the gas is very low (<1CV 7 ; 
e.g. ICaselli et all I1998D . shock waves develop a mag- 
netic precursor that accelerates, compresses and heats 
the ion fluid before the neutral one, forming a C- 
type shock fseelDraindfl980tlFlower fc Pineau des Foretsl 
Il995t iFlower et al.lll996t iLesaffre et al.ll2004ft . The sub- 
sequent velocity decoupling between the ion and neutral 
components of the gas produces the sputtering of the 
mantles of dust grains, injecting large amounts of molec- 
ular material into the far downstream gas (e.g. Draine, 
Roberge & Dalgarno 1983; Caselli, Hartquist & Havncs 
1997; Schilke et al. 1997). 

Since silicon is heavily depleted onto dust grains 
in dark clouds such as L183 and L1448 (SiO 
abundance of <10~^: [Ziurvs. Friberg fc Irvine! 119891 : 
iReauena- Torres et al.ll2007l) . the detection of SiO is a 
strong indicator of the dust grain processing by C- 
shocks. The SiO abundances in molecular outflows such 
as L1448-mm/IRS3 are >10~ 6 , and the typical SiO emis- 
sion arising from the warm postshock gas shows broad 
line profiles red- or blue-shifted with respect to the ve- 
locity of the ambient cloud (e.g. iMartin-Pintado et al.l 
fl99l . 
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iLefloch et all (|1998h and Uimenez-Serra et all (|2004l ) 
have reported the detection of very narrow SiO emission 
at ambient velocities with abundances of only ~lCP n - 
1(T 10 , toward the NGC1333 and L1448-mm molecular 
outflows. For L1448-mm, it has been proposed that this 
emission, which is also correlated with an enhancement 
of th e ion and electron densities ijJimenez-Serra et al.l 
arises from material recent ly processed by the 
magn etic precursor of C-shocks (|Jimenez-Serra et al.1 
l200l . 

Modelling of C-shock waves has traditionally focused 
on, and successfully reproduced, the SiO abundances 
and line profiles generated in the evolve d postshock gas 
(jSchilke et al.lll997tlGusdorf et al.ll2008[ ). but it has usu- 
ally ignored that molecular abundances and line profiles 
must change as C-shocks evolve in yery young outflows 
like th at powered by L1448-mm. iJimenez-Serra et al.l 
(|2008l hereafter J08) have shown that the observed trend 
for the enhancement of the SiO and CH3OH abun- 
dances from the prec ursor to the high velocity compo- 
nents in this outflow (| Jimenez-Serra et al.ll2005l ). can be 
explained by the progressive sputtering of the grain man- 
tles and cores as the shock evolves. However, no attempt 
has been made so far to predict the expected molecular 
line profiles, and their changes, produced by the propa- 
gation of shock waves through the unperturbed gas. 

Here we present the first results of the evolution of the 
line profiles of SiO, H 13 CO+ and HN 13 C, produced by 
the propagation of C-shocks. The derived line profiles 
strongly suggest that the narrow SiO lines and the ion 
enhancement are fingerprints of the magnetic precursor 
of shock waves. They also reproduce the broader moder- 
ate and high velocity emission components arising from 
the postshock gas in L1448-mm at later stages of the 
evolution of the shock. 

2. THE MODEL 
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Fig. 1. — Velocity of the ion and neutral fluids, v n and Vi, and 
temperature of the neutral gas, T n , for two representative C-shocks 
with i> s =30 and 70 kms -1 . The initial H2 density and magnetic 
field are no=10 5 cm -3 and Bo=450 fj,G. The filled dots and vertical 
dotted lines indicate the flow times for which we freeze the evolution 
of the line profiles in both shocks (Sec. 3 and Fig[2j. T n in our 
fit profile increases later in the shock compared to MHD models. 
Although this is not crucial for the magnetic precursor stage, it 
may significantly affect the derived SiO line profiles at stage 2 in 
the 70 kms -1 shock (see text in Sec. 3 for details). 



(|1997f ) and J08, the abundances of sputtered material 
from dust grains, and in particular of silicon, hardly 
change with the initial H2 density of the gas. On the 
contrary, these abundances are drastically enhanced for 
increasing shock velocities due to the larger ion-neutral 
drift velocities, \v n — Vi\, achieved in the shock. We have 
assumed that a small fraction of Si (~0.01% with respect 
to H2O; J08) is present in the grain mantles. Projection 
effects are also taken into account in the calculation of 
the SiO abundances by including the inclination angle 
of the outflow with respect to the line of sight, 6 (J08). 
We consi der that SiO is either dire ctly released from the 
mantles (|Jimenez-Serra et al.ll2005l ) or rapidly formed by 
the oxidation of S i sputtered from the cores (i<25yrs; 
iCabrit et al.ll200"7t ). The gas phase abundance of SiO, 
x(SiO), is calculated as a function of the flow time, 
t n , and the SiO column density in each plane-parallel 
slab of gas within the shock, iV(SiO), is estimated as 
7V(SiO)=x(SiO)n(H 2 ) Az, where n(H 2 ) and Az are the 
H2 density and spatial width of the slab, respectively. 

The level populations (i.e. excitation) of SiO in each 
slab of gas are calculated by using the LVG (Large Veloc- 
ity Gradient) approximation. We use the H 2 cq l lisiona l 
coefficients of SiO derived by iDavou fe Balancal (2006). 
The line profiles are then obtained by making the ra- 
diative transfer through all individual slabs of material. 
The line peak optical depth, tq, in each slab of gas is 
calculated by means of the LVG approximation through 
the following expression: 



We consider a plane-parallel C-shock propagating with 
a velocity v s through the ambient medium. The multi- 
layered steady-state physical structure of the shock in 
the frame of the preshock gas, is calculated by using 
the analytic approximation (fit) described in J08. Fig.Q] 
shows the shock field for v n , Uj and T n , as a func- 
tion of the flow time for two C-shocks with v s =30 and 
70 kms -1 , initial H 2 density n =10 5 cm -3 , and mag- 
netic field £? =450^G (B is calculated following the 

1 /2 

relation -60=^0( ^-3 ) with 6 ~1. See also Sec. 2 in 
J08 for the details on the approximation). The Alfven 
speed is 77^4=2. 18 kms -1 (derived from Bq and no as 
v A =—, — B ° ; see lMvers fc Goodmanl [l988). This ap- 

-v/47rm(//)no 

proximation has been validated in J08 (Sec. 4 in this 
work) and reproduces the physical structure of C-shocks 
for n =10 4 -10 6 cm . Hereafter, the term time evolution 
will be used for the variation of the molecular line emis- 
sion as a function of the flow time, which corresponds to 
the characteristic dynamical time across the shock. 

From Fig.[TJ it is clear that the ion fluid is accelerated 
before the neutrals by the magnetic precursor at the early 
stages of the shock. Since dust grains are charged par- 
ticles, they are coupled to the magnetic field and move 
with the ion velocity, Vi. In the postshock gas, the veloci- 
ties of the charged and neutral fluids re-couple again due 
to the transfer of momentum from the charged particles 
to the neutral fluid by collisions. 

The line profiles of species linked to the grain chemistry 
such as SiO, are derived by considering the sputtering of 
the grain mantles and cores by collisions between dust 
grains and neutral particles such as H 2 , He and other 
heavier species including CO. As shown bv lCaselli et al.l 



. = ( c \ 3 g up Am fNi N u \ 
° \Vui) 7T dv z /dz\gi g u J 

where c is the speed of light, and v u i and A u \ are the 
frequency and Einstein coefficient of the molecular tran- 
sition for which we calculate the line profiles. N u and Ni 
are the level populations of the upper and lower energy 
levels respectively, and g u and gi, the degeneracy asso- 
ciated with them. dv z /dz is the radial velocity gradient 
along the propagation of the shock. 

To determine the behaviour of the optical depth, t, 
as a function of velocity, we convolve the peak optical 
depth, ro, at a given velocity (slab of material) with a 
Gaussian-like profile. The linewidth of every individual 
profile, Av, is estimated by assuming thermal and tur- 
bulent broadening, Avth and Av tur - The thermal broad- 
ening depends on the temperature of the neutral fluid, 
T n , as: 

/ 8 ln2 k B T n 

Av th = \ (2) 

V m 

where ks is the Boltzmann constant and m is the mass 
of the molecule. The thermal broadening for the SiO 
emission at the precursor stage with T„=10K, is only of 
Au t/l =0.1kms -1 . 

The turbulent linewidth is considered as a fixed param- 
eter throughout the shock without largely affecting our 
estimates of the total linewidth. Althoug h turbulence 
in sho cks still remains poorly understood, iJoulan et al.l 
( 1998) studied the effects of viscous dissipation and ion- 
neutral velocity decoupling on the chemistry and molec- 
ular line emission arising from turbulent flows (filamen- 
tary vortex) in the presence of magnetic fields in diffuse 
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clouds. In their model, Uoulan et alj (|1998f ) determined 
that the maximum tangential velocity of the vortex fil- 
aments (which gives t he rms velocity dispersion of the 
large scale turbulence; I Jimenez] 1 19971 iBelin et aT1ll996D 
is close to the Alfven velocity in the neutral fluid, va- If 
we therefore assume that Ai> tur is directly proportional 
to va, the increase of Av tur in the postshock gas can 
be estimated from the compression of the magnetic field 
and of the density of the gas at the late stages of the 
shock. This compression is given by the shock jump con- 
ditions B^Bq(v s /vo) and n=no(v s /vo), where v s is the 
shock speed and vq depends on va and v s as described 
in J08. From here, Avtur in the postshock gas ranges 
from ^2.5 to 3.5 times larger than in the preshock gas 
for shock speeds between v s =30 and 70 km s _1 . However, 
despite the increase by a factor of ^3 in Avtur, the total 
line-width of the SiO line profile only differs by a factor of 
^1.2 from that calculated assuming constant turbulent 
broadening, because thermal broadening is dominant. 

For the ion and neutral species H 13 CO + and 
HN 13 C, we assume that their flux density, and there- 
fore, their abundances remain roughly constant across 
the shock. This is supp o rted by the results of 
iMarkwick. Millar fc Charrdevl (pOOOf ). in which HCO+ 
and HNC are chemically inert after the disruption of the 
grain mantles into the gas phase by the passage of a 
C-shock. The initial column densit ies for these molec- 
ular species are 7V~10 12 cm -2 ( Ji menez-Serra et al.l 
l200l . and as the line profiles of H 13 CO+ and HN 13 C 
get broader with the evolution of the shock, their 
peak column densities, Nq, are re-scaled following 
Nq—NSv* /Av*, where 6v* is the variation of the neu- 
tral/ion velocity from one slab of gas to the next one, 
and Av* is the total increase of the neutral/ion velocity 
from the beginning of the shock to the current slab of 
material in the shock. The excitation of H 13 CO + and 
HN 13 C has been derived by using the LV G approxima- 
tion a nd the H2 collisional coefficients of ISchoier et al.1 
(|2005l) . The line profiles are derived following the same 
procedure as that used for SiO. Finally, we assume that 
the emission of H 13 CO + and HN 13 C, as well as that of 
SiO, has the same spatial extent and morphology. Beam 
filling factors have not been considered in these calcula- 
tions. The line intensities therefore depend on the size 
of the shocked region. 

A bow-shock morphology for the shock may lead to 
substantial variations on the derived SiO line profiles 
with respect to those calculated in the plane-parallel 
case. The planar geometry is however easier to handle 
in the calculation, and since the SiO bullets in L1448- 
mm only show a partial bow-shock morphology (see 
iDutrev. Guilloteau fc Bachillerl [19971) . the use of both 
model geometries are equally justified. 

3. TIME EVOLUTION OF THE MOLECULAR LINE 
PROFILES IN C-SHOCKS 

In Fig.[2l we show the evolution of the line profiles 
of SiO, H 13 CO+ and HN 13 C obtained by freezing the 
propagation of two C-shocks, one with ^=30 km s _1 and 
the other with u s =70kms _1 , at 5 different flow times 
(see dots in Fig.[T]). In order to directly compare our 
results with observations, we have included the preshock 
gas velocity, v c i , in the calculation of the molecular line 
profiles. For the L1448-mm outflow (Sec. 4), we assume 



that i> e ;=4.7kms 1 an d the inclination angle is 0=70° 
(jGirart fc Acordl l20?)l . The SiO, H 13 CO+ and HN 13 C 
line profiles of Fig. [2] are then derived as a function of the 
radial velocity measured by the observer, vlsr- While 
for HN 13 C this velocity is Vlsr=v c i + v n , for H 13 CO+ 
v lsr is v c i + Vi. For SiO, one may think that this molec- 
ular species would move with the ion fluid once it is 
ejected from grains. However, since SiO is not coupled to 
the magnetic field, this molecule is rapidly slowed down 
by collisions with H2 to the velocities of the neutral fluid 
in time-scales, t s , of only few days. For a SiO molecule 

t s - 1 =n( J ff 2 ) 7 ra( J ff 2 ) 2 v /|^«2xl0- 12 T i y 2 nosec- 1 . 

Assuming n =10 5 cm- 3 , T n =10K and a( J ff 2 )=0.7414 A 
as the internucle a r dis tance for a H 2 molecule 
(|Huber fc Herzberel Il979f ). the slow down time is 
t s ^2Q days. 

From Fig. [21 we find that the SiO emission generated 
in shocks has very different line shapes depending on 
the evolutionary stage of the shock. In particular, it is 
striking that the typical SiO line profiles arising from 
the magnetic precursor at t„=3yrs, have very narrow 
linewidths of only ~0.5kms _1 and are slightly shifted 
in velocity (vlsr ^5.2 km s^ 1 ) with respect to the am- 
bient cloud (vlsr ~4.7kms~ 1 ). The narrow linewidths 
of this emission are due to the low temperatures of the 
neutral fluid at the magnetic precursor stage (T„=10K; 
see Fig.[T]), and the small velocity shift is explained by 
the fact that the sputtering of the grain mantles re- 
quires ion-neutral drift velocities of \y n — Vi |>6 km s _1 to 
generate SiO abundances of >10 -12 . This occurs when 
H 13 CO + reaches velocities of uls,r~12 kins -1 that im- 
ply \v n — «i|~7kms _1 (see Fig. [I]). The peak intensity 
of the narrow SiO line is 4 orders of magnitude larger 
for the ■u s =70kms- 1 shock (0.08K vs. 2xlO" 6 K) due 
to the larger SiO abundances (~10~ 8 vs. 10 -11 ) injected 
from the grain mantles into the gas phase. This intensity 
is increased by only a factor of 3 if we consider that the 
fraction of silicon within the mantles is of 0.1% with re- 
spect to H 2 (i.e. 1 order of magnitude larger than that 
assumed in Sec. 2). This is due to the fact that only a 
fraction of the grain mantles (<50%) is released into the 
gas phase for the same flow time of £„=3yrs. 

For later flow times, the narrow SiO feature is swept 
out by the larger efficiency of the sputtering of the grain 
cores once th e velocity threshold of >25km s -1 i s over- 
come (see e.g. ICaselli et al1fl997l lMa r TeTliI1 l2"000') . This 
enhances the gas phase abundances of SiO by several or- 
ders of magnitude. For the 30kms _1 shock, these abun- 
dances are of ^10 -8 and 10~ 7 at t n — 10 and 20yrs (see 
J08). For the 70kms~ 1 shock, the SiO abundances are 
enhanced up to 5x 10~ 6 and 7x 10~ 6 at t n =13 and 22 yrs, 
respectively. The increase of T n with the propagation of 
the shock leads to larger peak intensities and broader 
SiO line profiles. The SiO emission subsequently de- 
velops double peaked line profiles. The lower velocity 
peak is centred at vlsr. ~12kms _1 for v s =30kms -1 
and at 32kms~ 1 for w s =70kms~ 1 . This peak is cor- 
related with the maximum ejection of silicon from the 
grain cores (~10~ 7 and 7xl0~ 6 , i.e. 0.3% and 19% 
of all silicon locked in the cores for a SOkms" 1 - and a 
70kms _1 -shock, respectively), which happens as soon as 
the ion-neutral drift velocity reaches the threshold value 
of 25kms~ 1 . Since the fraction of silicon released into 
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SiO (2^1) H 13 C0 + (l->0) + HN 13 C (l->0) 

v a = 30 km s -1 v s = 70 km s _1 v s = 30 km s -1 v s = 70 km s" 




Fig. 2. — Time evolution of the SiO 2— >1, the H 13 CO + 1— >0 and the HN 13 C 1— >0 line profiles along the propagation of two C-shocks 
with t) s =30 and 70kms — 1 . Numbers in parenthesis indicate the stages at which the profiles are frozen (see Fig.[TJ. The flow times are 
shown in the upper part of all panels. Line intensities are in units of radiation temperature and the vertical dotted lines show the velocity 
of the preshock medium (4.7kms _1 ). 



the gas phase is larger for the 70kms _1 shock (by almost 
a factor of 100), the peak intensities therefore differ by 
more than 2 orders of magnitude from the i> s =30kms _1 
shock. We note that the mod e rate v elocit y peak was also 
presen t in the iSchilke et al.1 (|1997l ) and iGusdorf et al.1 
(|2008f) models, but averaged with and dominated by the 
brightest high velocity postshock emission. These mod- 
erate velo city profiles are sim ilar in shape to those re- 
ported bv lCabrit et all (120071) f or the young HH212 SiO 
microjet and by iLefloch et all (|1998fl in the NGC1333 
outflows. 

The high velocity peaks at v^sr ~30 and 70kms~ 1 
(i„^500yrs; Fig. [2]), are produced by the large accu- 
mulation of SiO in the postshock wide layers (slabs) 
of gas. Although the derived brightness temperatures 
are very high, this velocity component qu alitatively re- 
sembl es th e SiO line profiles ob tained by ISchilke et al.l 
(I1997D and Gusdor Tet alj (|2008l), and those observ ed by 
Uimenez-Serra et al.l (|2004l ) and lNisini et~aTl (|2007f ). The 
lack of SiO gas-phase chemistry or depletion onto grains 
in our model may account for the too strong emission 
predicted to arise from the high velocity postshock gas 
(see Sec. 4). 

We note that discrepancies in T n do exist be- 
tween the fit profile calculated with our approxima- 
tion and the solutions derived from comprehensive MHD 
models where T n increas es earlier in the shock (e.g. 
Kaufman fe Neufeldl H996L see also Fig. 2 in J08). This 
is not crucial for the precursor stage where MHD calcu- 
lations predict an initial lo w-temperature plateau for T n 
(with T„~10-20K; see e.g. iFlower fc Pineau des Foretsl 
|2003| ). However, for the moderate and high velocity com- 



ponents, the earlier increase of T n in the shock may pro- 
duce significant variations in the peak intensity of the 
SiO line profiles compared to our results. While for 
the SOkms^ 1 shock the SiO peak intensity could be in- 
creased by only a factor of 2, for the 70kms~ 1 shock 
this intensity could be enhanced by more than 1 order 
of magnitude. Therefore, although quantitative results 
cannot be given for the moderate and high velocity com- 
ponents, specially for the 70kms~ 1 shock (also note that 
out model lacks of SiO gas phase chemistry and freeze- 
out; see Sec. 4 for a more detailed discussion), the qual- 
itative behaviour of the SiO line profiles is preserved. 
Quantitative predictions of the SiO line fluxes across the 
shock will be carried out in the future with a compre- 
hensive MHD code including chemical processes. 

Regarding the broadening of the SiO lines, the earlier 
increase of T n in the shock does not substantially affect 
the derived linewidths of the SiO line emission. Indeed, 
for both cases with u s =30kms _1 and v s =70kms~ 1 , 
these linewidths differ by less than a factor of 1.4 from 
those shown in Figure[2] at every evolutionary stage 
within the shock. 

For H 13 CO + and HN 13 C, the differences in the line 
profiles of these species become apparent at the mag- 
netic precursor stage. The line profile of HN 13 C peaks at 
the preshock velocity of 4.7kms _1 at i„=3yrs (stage 1), 
and has a linewidth of ~0.7-0.8 km s -1 . The line emission 
of H 13 CO + , however, is (red)shifted to i>LSij=5.2kms -1 
and shows a moderate velocity wing with a terminal ve- 
locity (maximum velocity attained by the ion fluid at 
that certain flow time and reproduced in the molecu- 
lar line profile) of ^12kms~ 1 for w s =30kms _1 , and of 
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^23 kms -1 for i; s =70kms -1 (see Fig.[l|. We note that 
the shift in velocity of the H 13 CO + intensity peak is sen- 
sitive to the time step used in the calculation of the 
profiles. The line profiles of H 13 CO+ and HN 13 C get 
broader as the shock evolves. While the terminal veloc- 
ity of the H 13 CO + emission is reached at stage 2 (t n =2l 
and 13yrs; see Fig. [5]), the line profiles of HN 13 C require 
flow times of £„=500yrs to show broad emission with a 
similar terminal velocity. 

4. OBSERVATIONAL FINGERPRINTS OF THE MAGNETIC 
PRECURSOR OF C-SHOCKS 

In panels (a), (b) and (c) of Fig.|4j we show the 
comparison between the line emission of SiO, HN 13 C 
and H 13 CO+ observed toward the L1448 -mm outflow 
(|Jimenez-Serra et al.1 I2004T ) . and the line profiles ob- 
tained by our model for a C-shock with 7j s =70kms _1 , 
n =10 5 cm- 3 , £ =450 ^67, T =10K and v d =4.7kms -1 
at the magnetic precursor stage with £„=3yrs. The 
shock velocity of 70 kms -1 is selected based on the ter- 
minal velocity of the SiO emission in L1448-mm. Al- 
though the time-scales of 3 yrs at the precursor are very 
short, they are consistent with the dynamical ages of the 
very young shocks det ected in the vicinity o f the L1448- 
mm source (i~90 yrs; iGirart fc Acordll2001h. and in the 
HH212 SiO microjet f£<25 yrs: ICabrit et alJfeOOTf ). The 
freezing flow time of i n =3yrs implies a radial velocity 
of ULSfl=5-2kms~ 1 (i.e. v n —0.5 kms -1 ), at which the 
observed narrow SiO lines show their emission peaks to- 
ward L1448-mm. The predicted peak intensity for the 
J=2^1 transition is of 0.08 K and implies a SiO abun- 
dance of 2xl0 -8 . We note that we require a shock with 
i> s =70kms -1 in order to inject this SiO abundance at 
this early stage of the shock. One may think that SiO 
abundances of ~10 -8 could also be obtained at later 
flow times in the shock. However, the calculated line 
profile would be considerably broader and its emission 
peaks would be largely (red)shifted, which is not con- 
sistent with the narrow SiO emission observed toward 
this outflow. Therefore, the radial velocity and the nar- 
row linewidth of the narrow SiO emission are strong con- 
straints for the determination of the shock age associated 
with this emission. 

The derived peak intensity of the narrow SiO emission 
shown in Fig. [4] only differs by 20% from the observed 
one. The estimated J=3^2/J=2— >1 line intensity ratio 
(of ~0.4) is, however, a f actor of 3 smaller than that mea- 
sured in L1448-mm (see lJimenez-Serra et al.l l2004). sug- 
gesting that the H2 densities toward this outflow might 
be larger than those considered in our model. The low ve- 
locity resolu tion (of >1 kms" 1 vs. 0.14 kms -1 shown in 
Fig.rjJ of the lNisini et al.l (|2007f ) observations toward this 
outflow, prevents the direct comparison of these single- 
dish data with our results for high- J SiO transitions. In 
any case, the derived line intensities of the SiO transi- 
tions with J up >5 are expected to be far below the 3cr 
noise level of these observations. 

The SiO gas phase abundance required to match 
the intensity of all narrow SiO lines is of ^2xl0 -8 , 
three orders of magnitude larger than that derived by 
iJimenez-Serra et al.l (|2004| ) . Beam dilution could explain 
the discrepancies in the SiO abundances. If the emitting 
region of the narrow SiO J=2^1 emission is <1" (of 
^1.5-3x 10 15 cm for the magnetic precursor at the dis- 



tance of 200 pc; see J08), the predicted abundance of 
~2xl0 -8 is diluted to ~10 -11 , which is of the same 
order of magnitude as th at measured in this outflow 
pimenez-Serra et al.ll2004f ). 

From Fig.lU we also find that the difference between 
the H 13 CO + and the HN 1 3 C line profiles observed by 
IJimenez-Serra et al.l (|2004f) toward L1448-mm, are re- 
produced by the model. While the neutrals (HN 13 C; 
panel b) remain unperturbed at the early stages of the 
shock (their line profiles are centred at the velocity of 
the preshock gas of 4.7 kms -1 ), the ion H 13 CO + fluid 
(panel c) has slipped to larger velocities (with its line 
peak centred at ~5.2kms -1 ) and shows an ion enhance- 
ment in the precursor compo nent that is consistent wit h 
that measured in L1448-mm (|Jimenez- Serra et al.ll200l . 

In contrast with the observed line profile of H 13 CO + , 
however, the calculated line of this species shows 
redshifted broad emission with terminal velocities of 
<12kms -1 (red cl histogram; Fig.|4j. The obser- 
vational procedure, frequency-switched mode, could 
be the responsible of the missing broad emission in 
the final spectrum of H 13 CO + . Another possibil- 
ity could be that H 13 CO + rapidly recombines in the 
magnetic precursor due to the large electron density 
enhanc ement at this stage of the evo lution of the 
shock (jFlower fc Pineau des Foretsl 12001 . If we as- 
sume that the H 13 CO + abundance decreases by ion re- 
combination with the flow time (we consider a frac- 
tional i onisation of ~5xl0~ 5 in the precursor com- 
ponent; IJimenez-Serra et al.l l2006h . the broad H 13 CO+ 
wing emission vanishes, leading to a narrower H 13 CO + 
line profile that mimics the observed one (terminal veloc- 
ity of ^7 kms -1 ; blue c2 histogram). In the final derived 
spectrum of H 13 CO + , we have also included the ambi- 
ent cloud component at 4.7 kms -1 with a H 13 CO + peak 
colum n density of 5xlO n cm -2 (|Jimenez-Serra et al.l 
120041 ). The calculated line intensities of H 13 CO+ and 
HN 13 C differ by 13-50% from those measured toward 
L1448-mm. Variations by a factor of 2 in the molecu- 
lar column densities and initial H2 density would lead 
to H 13 CO+ and HN 13 C line intensities that differ by 
less than a factor of 2 from the ones calculated for 
iV(H 13 CO+)=A(HN 13 C)=10 12 cm -2 and rj =10 5 cm -3 . 
Interferometric observations that provide the morphol- 
ogy and spatial extent of this emission are therefore re- 
quired to solve the discrepancies in both the molecu- 
lar abundances and line peak intensities of H 13 CO + and 
H N 13 C in the magneti c precu rsor. 

IJimenez-Serra et all ((2004) discussed the possibility 
that narrow SiO could be generated by the evap- 
oration of the icy grain mantles by the illumina- 
tion of U V radiation in Photon Dominated Regions , 



orPDRs (Walmsle y. Pineau des Forets fc Flowed 1999) 



iSchilke et all (|2001h reported the detection of narrow 
SiO emission with SiO abundances of ~10 -11 toward a 
sample of PDRs such as the Orion Bar or S 140. The 
linewidth of this emission is however of ^2-3 kms -1 , a 
factor of ^5 broader than those measured toward L1448- 
mm. Furthermore, HCO, one of the main products of 
the UV photo-che mistry (see e.g. lSchenewerk et al.lll988l : 
iGerin et"aT]l2008h is not detected in the precursor com- 
ponent of this outflow, which allows to rule out the PDR 
origin for the narrow SiO emission. 
Alternatively, narrow SiO emission may arise from an 
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Fig. 3. — Comparison between the predicted line profiles of SiO, HN 13 C and H 13 CO + at t n =3yrs for a shock with u a =70kms~ 1 and 
rao=10 5 cm - 3 (upper panels), and the line profiles observed toward the L1448-mm outflow (lower panels; Jimcncz-Scrra ct al. 2004). The 
ambient (4.7kms _1 ) and the precursor (5.2kms _1 ) components are indicated by the vertical dotted lines. In the upper (c) panel, we show 
the line profiles of H 13 CO + calculated with and without ion recombination (blue c2 and red cl lines, respectively). 



- 0.05 




Fig. 4. — SiO 2— >1 line profiles (black solid lines) derived for 
the same C-shock model (t) s =70kms — 1 and no=10 5 cm -3 ) shown 
in Figs. [2] and [4] which best fit the full observed SiO line profile 
toward L1448-mm (lower panel; Jimcncz-Scrra ct al. 2004). The 
freezing flow times (t n =3, 11 and 440 yrs) are shown in every panel. 
The vertical dotted lines indicate the radial velocities (i;^gfl=5.2, 
12 and 68kms — 1 ) that correspond to these flow times. We also 
include the results of a second C-shock model that considers the 
possible effects that gas-phase chemistry may have on the derived 
SiO line profiles (red thick solid lines). 



intervening cloud where turbulence motions of at least 
6kms _1 would be efficient enough to erode the grain 
mantles and inject detectable abundances of SiO into 
the gas phase. However, in contrast with the large- 
scale morphology expected from this scenario, the nar- 



row SiO emission is not widespread around L1448-mm 
as shown by the low upper limits of the SiO abundance 
of <10~ 12 found in the regions wit h no p recursor in- 
teraction ()Jimenez-Serra et all 12004 120051) . We note 
that these upper limits are of the same order of magni- 
tude as those measured in the quiescent gas of the L183 
dark cloud by the high-sens itivity 30 m observations of 
iReauena- Torres et al.l (|2007|) . This indicates that the 
lack of SiO emission toward these regions is caused by 
a real lack of SiO in gas phase. Prom all this, we pro- 
pose that the most likely explanation for the narrow SiO 
lines detected toward L1448-mm, is the interaction of the 
magnetic precursor of C-shocks in very young outflows. 

The line profiles calculated across our parametric C- 
shock structure also suggest that the moderate and high 
velocity SiO components observed toward L1448-mm are 
linked to later stages in the evolution of the shock. In 
Fig.QJ we show the comparison between the full SiO 
line profile in L1448-mm (lower panel), and the SiO 
line emission predicted for a shock with w s =70kms _1 
and no=10 5 cm -3 at three different evolutionary stages 
(black solid lines). The narrow SiO emission is exactly 
the same as that shown in panel (a) of Fig.[4]for i n =3 yrs, 
and the moderate velocity component is reproduced by 
a C-shock at t n =llyrs. We note that the discrete na- 
ture of these two components is a direct consequence 
of the assumed chemical structure of dust grains: icy 
mantles surrounding the silicate/graphite cores, which 
have different binding energies (of 0.53 eV and 50 eV for 
the mantles and the co res, respectively; see lTielens et al.l 
11994 iMav et all 12000ft . The high velocity gas emission 
is matched with a C-shock at i„=440yrs. 

Although the model fails to reproduce the peak inten- 
sity of the moderate and high velocity SiO components, 
the line shapes are qualitatively similar to those observed 
in L1448-mm. At i„ =440 yrs, however, the moderate ve- 
locity emission shows an excess in intensity which is dif- 
ficult to match with the observations. The facts that de- 
pletion of Si/SiO onto dust grains and/or that silicon gas- 
phase chemistry have not been taken into account in our 
model, could ex plain these discrepanc i es. In cluding the 
depletion rate of lTielens &: H ollcnbacbl f| 19851 . see Eq. 4 in 
this work) with a sticking coefficient of ~1, the depletion 
of SiO onto dust grains reduces the abundance of this 
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molecule by a factor of *~1.2 for intermediate flow times 
and by a factor of ~5 in the far postshock gas. More 
dramatically, the conversion of silicon into SiO via gas- 
phase reactions with OH and O2 (e.g. Si+OH — » SiO+H 
and Si+02 — > SiO+O) is known to be not 100% efficient 
for high densities (of no=10 5 -10 6 cm -3 ), for which OH 
and O2 are partially de stroyed by the atom ic hydrogen 
generated in the shock ([Gusdorf et al.l [20081 ). This may 
lead to significantly lower SiO abundances (by more than 
1 order of magnitude; see lGusdorf et al.l|2008f ) than those 
obtained directly from the sputtering of dust grains. 

If we now consider that the gas-phase abundance 
of SiO cannot exceed ^10~ 6 (note that this is the 
maximum SiO abundance estimated in L1448-mm; see 
iMartm-Pintado et al.lll992l ). the SiO line profile derived 
at £„=440yrs mimics the emission observed toward this 
outflow (see red solid line in panel (c) of Fig. 2]). The 
decrease of the gas-phase abundance of SiO does not 
largely affect the SiO lines calculated for the precursor 
and moderate velocity components at t n =3 and llyrs, 
because the grain core sputtering is not efficient enough 
to inject SiO abundances larger than 10~ 6 (see panels 
(a) and (b) in Fig.[4j. Beam filling factor corrections 
could also account for the different peak intensities ob- 
tained from the model for the moderate and high ve- 
locity emission of SiO, assuming that both components 
show different spatial extents. High-angular resolution 
observations with interferometers such as the PdBI and 
the SMA, are therefore strongly required not only to bet- 
ter constrain the initial conditions and the evolutionary 
stage of the different shocks suggested by the model, but 
also the chemical properties of the dust grains present in 
this star forming region. 

Finally, shock velocities of 70kms~ 1 clearly exceeds 
the critical velocities of C-shocks (of ^45-50 km s -1 ; 



iDraine et all Il983f ). iLe Bourlot et all (|2002f) have 
shown that these critical velocities can be increased to 
few ^lOOkms -1 for moderate H2 densities and high 
magnetic fields (of few mG) as those observed to- 
ward low mass star forming regions such as NGC1333 
(|Girart et alj 120061 ). However, we cannot rule out 
the possibility that the high velocity SiO emission 
toward L1448-mm is associated with a J-type dis- 
continuity within the postshock gas in this outflow 
(jChieze. Pineau des Forets fc Flowerlll998f) . 

In summary, we have modelled for the first time the 
expected evolution of the line profiles of SiO, H 13 CO + 
and HN 13 C as C-shocks propagate through the ambient 
medium. The model suggests that the narrow SiO emis- 
sion and the ion enhancement found in the young shocks 
of the L1448-mm outflow, are signatures of the magnetic 
precursor of C-shocks. The different velocity components 
(precursor, moderate velocity and high velocity gas) ob- 
served toward this outflow could be explained by the co- 
existence of C-shocks with w s =70kms _1 at different evo- 
lutionary stages within the single-dish beam of the SiO 
observations. The model will be used in the future to 
reproduce the molecular line emission observed toward 
a sample of molecular outflows under different physical 
conditions and different evolutionary stages. 
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